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Abstract The Lunar CRater Observations and Sensing Satellite (LCROSS) mission 
impacted a spent Centaur rocket stage into a permanently shadowed region near the 
lunar south pole. The Sheperding Spacecraft (SSC) separated ^9 hours before impact 
and performed a small braking maneuver in order to observe the Centaur impact plume, 
looking for evidence of water and other volatiles, before impacting itself. 

This paper describes the registration of imagery of the LCROSS impact region from 
the mid- and near-infrared cameras onboard the SSC, as well as from the Goldstone 
radar. We compare the Centaur impact features, positively identified in the first two, 
and with a consistent feature in the third, which are interpreted as a 20 m diameter 
crater surrounded by a 160 m diameter ejecta region. The images are registered to 
Lunar Reconnaisance Orbiter (LRO) topographical data which allows determination 
of the impact location. This location is compared with the impact location derived 
from ground-based tracking and propagation of the spacecraft's trajectory and with 
locations derived from two hybrid imagery /trajectory methods. The four methods give 
a weighted average Centaur impact location of -84.6796°, -48.7093°, with a la un- 
certainty of 115 m along latitude, and 44 m along longitude, just 146 m from the 
target impact site. Meanwhile, the trajectory-derived SSC impact location is -84.719°, 
-49.61 °, with a la uncertainty of 3 m along the Earth vector and 75 m orthogonal to 
that, 766 m from the target location and 2.803 km south-west of the Centaur impact. 

We also detail the Centaur impact angle and SSC instrument pointing errors. Six 
high-level LCROSS mission requirements are shown to be met by wide margins. We 
hope that these results facilitate further analyses of the LCROSS experiment data and 
follow-up observations of the impact region. 
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1 Introduction 



On October 9, 2009 at 11:31:19.506 UTC, a spent rocket upper stage, the Centaur, 
impacted a Permanently Shadowed Region (PSR) within the crater Cabeus near the 
lunar south pole. The impact was observed by the Lunar CRater Observations and 
Sensing Satellite (LCROSS) mission's Shepherding Spacecraft (SSC) which carried 
nine instruments, including cameras, spectrometers and a radiometer. The impact was 
also observed from the Earth and by the Lunar Reconnaissance Orbiter (LRO) mission. 
The LCROSS experiment's primary objective was to determine the presence or absence 
of water ice in a permanently shadowed crater at the Moon's South Pole and to help 
constrain competing water distribution models. 

After launch, the LCROSS SSC controlled the Centaur for the next four months, 
performing maneuvers to target the Centaur at the planned impact site. Approximately 
nine hours prior to impact, the SSC separated from the Centaur, performed a braking 
burn to build a four-minute separation between itself and the Centaur at impact, and 
oriented its instruments toward the planned impact site within Cabeus. Fifty minutes 
prior to impact, the nine LCROSS instruments were powered on and began taking 
data. 

Initial findings based on the LCROSS data relating to the water measurements can 



be found in A. Colaprete et al. (2010), the impact cratering process is discussed in 



Schultz et al. (2010), and related results from LRO are in accompanying papers, 
e instruments are described in more detail in |Ennico, K. et al. ( 2008| 2010a|b|); pT] 
Idmann et al. (2010). Flight operations are described in |Tompkins, P. et al. (|2010a|b ); 
Strong, J. et al. |(2010), the trajectory and mission design in Cooley et al.| (2010). 
Overall program mission lessons learnt can be foun d in Andrews (|2010|). Additional 



detail appears in the NASA Flight Director's Blog ( [LCROSS Flight Director's Blog 
[2009 ). 

In this paper, we describe a set of interrelated post-mission analyses that provide 
context for the instrument data, focusing on registering imagery of the impact region 
and identifying the crater and its features. The next section summarizes the results 
from these analyses, subsequent sections then providing the full analysis. 



2 Summary of Results 

Figure 1 shows the impact crater in Mid-Infrared (Thermal) Camera #1 (MIRl), Near- 
Infrared Camera #2 (NIR2) and a ground-based X-band (3 cm) radar image from 
Goldstone Solar System Radar (GSSR) (which has been registered to the LRO topo- 
graphical data). The relative alignments of the MIRl and NIR2 images, taken shortly 
before the SSC impact, are well known due to ground calibration and registration of 
earlier imagery taken during the mission. The MIRl shows a hot spot approximately 
90 ± 25 m wide. The NIR2 image shows a feature in the same location, centred ap- 
proximately —84.6774, —41.691 degrees in Mean Earth, selenocentric coordinates. It is 
a a 62 ± 20m diameter dark region, surrounded by a lighter ring that is 158 ± 40m 
wide. This is interpreted to be the crater and the eject a blanket. Note that the esti- 
mated Centaur crater size is 25-30m based on impact modelling and constrained with 
LCROSS imagery dP. Schultz et al.|2010[ ). The ring is brightest in the East (and slightly 
south). The directional bias of the eject a ring and the impact angle with respect to the 
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Fig. 1 Left to right: LCROSS-MIRl, LCROSS-NIR2 and Goldstone radar images of the 
Centaur impact crater. These are co-registered (see Section 3) and to the same scale. The 
yellow circles mark regions centred on the impact location as estimated by the trajectory 
information alone (see Section 4), with a radius equal to the 3a errors of that method (423m). 
The black dot marks its centre. 

local terrain are both predominantly in the easterly direction (see Section 7) suggesting 
one possible explanation that the impact angle caused the eject a bias. 

A radar feature is seen in the GSSR also - this is an area of increased radar 
brightness around 30% smaller than the NIR2 feature, is thought to be due to surface 
roughness associated with the ejecta blanket surrounding the crater created by the 
impact, since it is much larger than the crater itself. Table [l] summarizes the averaged 
location and basic features for the Centaur and SSC impact sites and other results 
from this paper. Subsequent sections detail the derivation of each quantity. 



Feature Characteristic 


Centaur Impact 


SSC Impact 


Time (UTC) 


11:31:19.506 


11:35:36.116 


Lattitude (°) 


-84.6796 ± 0.004 


-84.729 ±0.001 


Longitude ( °) 


-48.7093 ±0.016 


-49.61 ±0.014 


Altitude (m) 


-3, 825 ±10 


-3, 809 ± 10 


Error in Lat. x Long, (m) 


119 X 46 


31 X 39 


Dark Region Diameter (m) 


62 ±20 




Light Ring Diameter (m) 


158 ± 40 




Impact Angle ( ° ) 


3.67 ±2.3 




Distance to Target (m) 


146 


766 



Table 1 Key results of locations estimates and impact craters for the SSC and Centaur. Errors 
given are Icr. Altitude is with respect to a mean lunar radius of 1737.4 km, and latitude and 
longitude are given in Mean Earth, selenocentric coordinates times are in UCT on 9th October 
2009. 



3 Registration of Key Images and Impact Feature Identification 

This section describes the process of registering images of the Centaur impact region 
against each other and the identification within the imagery of the impact feature. 
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Passive imaging of the impact site is difficult because it lies in permanent Sun-shadow. 
However, the MIR cameras could see the impact location due to the thermal signature 
of impact heated regolith. In addition, sufficient detail for registration is visible in 
some NIR2 camera images shortly before the SSC's impact because at this point the 
camera's field of view was contained fully within the PSR so the dynamic range of each 
image is reduced and its signal-to- noise ratio is improved. 

Naturally, active imaging methods do not suffer from the same issue and thus LI- 
DAR images from LRO (from the mini RF instrument) and GSSR images are available. 
The latter is possible because although in permanent Sun shadow, the floor of Cabeus 
crater is occasionally in Earth view. In total four data sets are used in the registra- 
tion here: imagery from LCROSS's Mid-Infrared Camera 1 (MIRl) and Near-Infrared 
Camera 2 (NIR2) shortly before SSC impact; Goldstone radar before (May 2009) and 
after (November 2009) impact observations; and LRO's Lunar Orbiting Laser Altime- 
ter (LOLA) topographical data. (Note that LRO mini-RF data were also registered but 
due to its lower spatial resolution compared to GSSR, it is not shown in the present 
analysis.) 

The registration process follows five steps. The first step is to identify the hotspot 
in the MIRl image. This is relatively simple because the average temperature of the 
surface within the PSR is considerably below the camera's detection threshold, residual 
heat from the Centaur impact is the only visible feature that moves from image to image 
- as the spacecraft frame moves over the lunar surface - over the last 30 seconds before 
communication was lost with the SSC. All other features remain fixed, e.g., hot pixels 
and faint ghosting from earlier images that persists until the next, automated flat field 
correction within the camera. 

The second step is to identify the Centaur crater feature in the NIR2 image. The 
relative pointing of the NIR2 and MIRl cameras are known from pre-launch calibration 
and, in addition, can be verified and updated by registration of overlapping images post 
launch (the latter method was used as the basis of alignment herein). However, because 
the LCROSS camera exposures were not time synchronized, we take the MIRl images 
immediately before and after each NIR2 image, identify the hotspot s in each, and 
interpolate between them to create a hotspot location corresponding to the time of the 
NIR2 image. In the case of the NIR2 image shown in Figure 1 (timestamp 3460653481, 
corresponding to 11:35:31.401 on 09 Oct, 2010, approximately 4s prior to it's own 
impact), the interpolated hot spot in the MIRl image lies directly on top of the impact 
feature seen in the NIR2 image, as is clear in the Figure. 

The third step is to register the NIR2 image to the GSSR radar data (in particular 
we use the exposure taken before the Centaur impact). To do this, the NIR2 image 
was approximately aligned with the GSSR image based on the spacecraft position and 
attitude that the time of the image was taken, together with camera pointing on the 
spacecraft, providing approximate latitude and longitude locations of the corner pixels. 
Then, eight features were identified that were present in both the NIR2 and GSSR 
images. Markers were placed in the centre of these features, and an elastic warping 
process (size, orientation, stretch and deformation) was then used to shape the NIR2 
image onto the GSSR image. This was done with the BunwarpJ plug in for ImageJ 
( I. Arganda-Carreras et aL]|2006 ). After the images were registered, more than seven 
additional features were found to be common to both images, giving confidence in 
the resulting registration. This was repeated now including all 15 common features as 
markers allowing a more precise registration. Again a further 5 features were found to 
be in common. Figure 2 shows the registered images and aligned features. 
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Fig. 2 GSSR X-band (3 cm) radar imagery overlaid on LRO-LOLA map (left) and with both 
the GSSR and LCROSS-NIR2 image overlaid (right). The blue arrows represent markers used 
for registration. The green arrows are additional features found to line up upon registration. 
The yellow circle represents the 3a error around the location estimate of the Centaur impact 
(upper right) and SSC impact (lower left), each based on trajectory data from tracking (Section 
4). These errors correspond to a radius of 423m in both latitude and longitude for the Centaur 
and 9x225m in the Earth vector (projected onto the lunar surface at the impact location) and 
orthogonal to the Earth vector for the SSC. The grid lines are at 1 ° intervals in longitude 
3.03 km) and 0.1 ° intervals in latitude (~ 2.82 km). 

The fourth step is to place the combined registered image set (MIRl, NIR2 and 
GSSR imagery) onto the LOLA terrain. This is straightforward since the GSSR data 
has already been registered to the LOLA data, and corner pixel coordinates were used 
to align the image. 

The fifth step is to identify the impact feature in the data sets. To do this, the 
trajectory-based impact location (see Section 4) was marked onto the registered set of 
images. An ellipse was drawn to mark the 3-sigma errors associated with the trajectory 
estimate. The MIRl image clearly shows a hot spot near the centre of the trajectory 
ellipse (see Figure 1). Moreover this is the aforementioned one feature that moves from 
image frame to frame, thus confirming it is the heat signature of the Centaur impact 
rather than noise. 

Precisely aligned with this in the NIR2 image is a feature described in Section 
2. For the detection in the GSSR data, the process was slightly more involved. The 
data were processed through a backprojection image formation processor to 5 m range 
resolution and 10m cross-range or azimuth resolution. The data were orthorectified 
during image formation by backprojecting onto the digital elevation model of the moon 
from the Japan Space Agency LALT laser altimeter instrument on the Kaguya lunar 
mission. The data were formed on an orthographic projection grid (note: this is not 
polar stereographic) . A region 15 km on a side was extracted from the pre-impact (May 
2009) and post-impact (November 2009) data, approximately centred on the LCROSS 
impact site. These single-look complex (SLC) data were averaged incoherently in range 
by a factor of 2 to form square (10 m by 10 m) pixels. These imagery were further de- 
speckled by convolving a unity gain boxcar filter 40m on a side (or 4 by 4 pixels at 
10 m resolution) with the data. The resulting imagery, now at approximately 40 m 
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Fig. 3 Goldstone X-band radar image of the Centaur impact region taken post impact and 
backprojection processed as described in the text. The blue spot shows the trajectory impact 
location as determined by a CFAR detection algorithm described in detail in the text. The red 
pixels represent areas which the CFAR program detected to be above its threshold in the ratio 
change image (constructed from before and after images). The red pixels are predominantly 
associated with the di?erence in shadows due to the di?erent viewing angles for the before 
and after impact exposures, however there is a small feature approximately 75m from the 
trajectory-based impact location. This feature has a 83% condence level and has a location 
consistent with the Centaur impact. 



resolution (but still at 10 m posting) were input to a constant false alarm rate (CFAR) 
thresholding program. A ratio change image was formed by dividing the November 
data by the May data, pixel by pixel. The resulting ratio image was fit to a log-normal 
probability distribution. The region chosen for this fit was the lower half of the image, 
or a 7.5 by 15 km region, which did not contain the putative impact zone. Once the 
log-normal parameters were extracted, an operating threshold of 10 false alarms per 
square kilometer was chosen, and the ratio image thresholded to highlight those pixels 
which contained statistically significant differences in reflectivity. This CFAR detector 
highlighted a point 75 m away from the impact point indicated by radio tracking of 
the final trajectory of the LCROSS impactor. From 2-D probabilistic analysis, for 10 
FA/km2, the probability that a random event lands within 75m of the impact zone is 
17%. Restating this result, there is 83% confidence in the result, or approximately a 
1 out of 5 chance that the detected point was a random event. The processed post- 
impact GSSR data around the LCROSS impact region is shown in Figure 3 (note: this 
image is different from the GSSR data shown in the other figures since it had different 
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-84.5 



-84.6 



-84.7 



-84.8 




Fig. 4 LCROSS-MIRl, LCROSS-NIR2, Goldstone Radar and LRO-LOLA images of the im- 
pact site co-registered. The grid hnes are at 1 ° intervals in longitude (~ 3.03 km) and 0.1° 
intervals in latitude (~ 2.82 km). The sun and Earth vectors at the time of the Centaur impact 
are shown (note that these vectors are different for the radar image which was taken about a 
month after the impact. The yellow circles and the arrows are the same as for Figure 2. 



processing, as described above). This cannot be taken as a statistically significant 
confirmation that this is the Centaur crater but it is consistent with it. 

Further analysis of the GSSR data could lead to a more denitive assessment of this 
feature. Indeed, to date only one exposure before and one after have been processed, 
whereas six before and four after exposures were taken. Combining these could improve 
the signal-to-noise ratio (particularly the radar speckle component) and allow clearer 
identification of the change feature. Further, of the available imagery, only the radar 
can potentially contain an SSC crater feature and thus this may enable verification of 
the SSC location. 

The complete set of MIRl, NIR2 and GSSR images registered to the LOLA back- 
ground is shown in Figure 4. Four factors increase our confidence in the image alignment 
and the location of the Centaur impact feature: (1) the alignment of ^ 20 features be- 
tween any two registered image pair (NIR2-GSSR, NIR2-L0LA and GSSR-LOLA); (2) 
the good alignment of the Centaur impact feature shown in the NIR2 and the MIRl 
hot spot; and (3) the radar feature in the 3a error margins of the trajectory-derived 
impact location of the Centaur (albeit not statistically significant). The agreement of 
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Table 2 Methods for locating the Centaur impact site 
Name Description 

Trajectory Only Based solely on orbit determination from the Deep Space Network 

(DSN) tracking and propagation in the lunar gravitational field 
Imagery Only Registration of NIR2 and GSSR imagery that has the impact 

Crater in view against LOLA-derived reference imagery 
Hotspot Forward Hybrid approach (manual): hotspot forward projection to the lu- 
Projection nar surface of the impact location hotspot in MIRl as determined 

by SSC position, attitude and camera pointing data 
Bundle- Hybrid approach (automated): bundle- adjusted (BA) image regis- 

Adjustment t rat ion against a LOLA-derived reference prior to impact, utilizing 

the trajectory and spacecraft attitude data 



this combination of imaging types increases our confidence in the identification of the 
impact site. 



4 Centaur Impact Location 

The Centaur impact location was estimated using the four, partially-independent tech- 
niques summarised in Table [2] Here we introduce all four methods and their results. 
A complete list of the data used and an overview of the information flow for each of 
these methods is given in Appendix A. Each method of impact location determination 
is described in full detail in Appendix B. In summary, the image only approach comes 
directly from the impact registration of Section [S] that provides a location of features 
seen in the images on the lunar surface. The trajectory only method uses just the DSN 
tracking data (Doppler and range data) along with propagation of that trajectory in 
the lunar gravity field, and other necessary perturbing forces to predict an impact 
location on the lunar surface, and no imagery at all. The hybrid approaches use the 
same imagery but not the registration itself, together with the trajectory data in or- 
der to estimate the impact location. Appendix B includes full details of the trajectory 
states (and propagation techniques) used that may be of use to other research efforts, 
whilst further information on the trajectory design and optimisation techniques use in 



LCROSS can be found in [Cooley etal] ( [20Tq| . 



Table [3] summarizes the results from each of the four methods, including the la 
errors. The final weighted average location is estimated to be -84.6796 degrees north, - 
48.7093 degrees west, in Mean Earth, selenocentric coordinates. This equates to 106.074 
km, -120.85 km, -1726.1 km in Cartesian coordinates. Figure 5 shows these results for 



Table 3 Summary of the Centaur impact location, using four different methods. Altitude is 
with respect to a mean lunar radius of 1737.4 km, and latitude and longitude are given in 
Mean Earth, selenocentric coordinates. The trajectory estimate given is that which includes 
the outgassing adjustment (see Appendix B) 



Method 


Lat (°) 


Long (°) 


Alt. (km) 


Error (m) 


1. Image Only 


-84.6774 


-48.691 


n/a 


130x130 


2. Trajectory 


-84.6827 


-48.6688 


-3.829 


141 X 141 


3. Hybrid: Manual 


-84.6805 


-48.7399 


n/a 


45 X 87 


4. Hybrid: Auto. 


-84.6782 


-48.7214 


-3.816 


66 X 102 


Weighted Ave. 


-84.6796 ± 0.004 


-48.7093 ± 0.016 


-3.825 ±0.01 


115 X 44 


Target 


-84.675 


-48.725 


-3.827 
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Easterly Displacement from Target location |km| 



Fig. 5 SSC and Centaur impact locations estimates compared. The two figures show the 
clusters of estimates for the SSC (left) and Centaur (right) using their 2 and 4 methods 
respectively, as well as the original target coordinates. All distances are in kilometers. The 
origin of both images is the Centaur impact target location. The ellipses represent la errors 
for each method - each centred on their respective method, apart from the Centaur trajectory 
method without outgassing whose error ellipse (pale blue) has a bias towards the south-south- 
east (the anti-sunward direction). The right hand image also shows the overall error ellipse in 
green centred on the weighted average impact location. 

the Centaur on the right and, separately, for the SSC on the left. (The process for the 
SSC is described in Section 4). 

The distance from the weighted average location to the target location is 139 m 
North and 44 m East, 146 m in total. The la error on this is 115m along latitude and 
44m along longitude. Thus the worst case mean plus 3a range is < 500 m. 



4.1 Centaur Location Error Analysis 

The four methods were combined by a weighted average to generate the impact location 
in Table [3] above. The relative weights are inversely proportional to the error of each 
measurement. The errors given are la assuming a normal distribution of measurements, 
and equate to a distance of 115 m from the mean position in the North-South direction 
and 44m in East- West. An alternative statement of the errors is that all estimates lie 
within a 148m radius circle centred on the average. 

Figure 5 makes clear that all the location methods are just about mutually consis- 
tent at the la level. A tightening of the error ellipsoids may be possible from a variety 
of sources: (1) imagery of the crater from LRO/LROC, LOLA or MiniRF; (2) improved 
imagery from GSSR allowing confirmation of the Centaur (and even SSC) impact fea- 
tures in that data; (3) image bundle adjustment registration which also refines the 
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Table 4 Summary of the SSC impact location, using two different methods. Latitude and 
longitude are given in Mean Earth, selenocentric coordinates 



Method 


Lat (°) 


Long (°) 


Alt. (km) 


Error (m) 


Trajectory 


-84.719 


-49.61 


n/a 


3 X 75 


Target 


-84.729 


-49.36 


-3.8091 





spacecraft trajectory and attitude and thus updates the trajectory based approaches; 
or (4) updated timing information on the spacecraft electronics boards. 



5 SSC Impact Location 

The SSC impact location was only possible to estimate using the trajectory method. 
The final weighted average SSC location is estimated to be -84.719 degrees north, - 
49.61 degrees west, in Mean Earth, selenocentric coordinates. This equates to 103.364 
km, -121.54 km, -1726.21 km in Cartesian coordinates. The results are summarised in 
Table [4] The analysis is described in full detail in Appendix C. 

The SSC was tracked until 1 hour before its impact, whereas the Centaur was 
untracked after separating from the SSC, thus the SSC has more precise trajectory data. 
However, the Centaur crater was imaged by the SSC cameras and as a consequence 
additional methods could be used to locate the Centaur impact, while only one was 
used in to locate the SSC impact (see Appendix A). Also, the Centaur crater produced 
a larger crater, easier to see in the radar image. 

The distance from the trajectory location to the target location is ^^766 m. The 
errors on the SSC trajectory as calculated by the JPL Orbital Determination team 
are 3 m x 75 m, Icr, where the 3m is in the Earth vector (as projected onto the lunar 
surface) and the 75 m is orthogonal to that. These are dominated by the errors of the 
tracking process. Thus the worst case mean plus 3a range from target impact location 
is ~ 1000 m. This is 2.803 km south-west of the Centaur impact. 



5.1 SSC Impact Error Analysis 

The averaged SSC location is almost 800 m from the target SSC location. This is 
compared to 5x smaller distance for the Centaur estimated impact point from its 
target. This can be explained by the fact that during the final mission events, and in 
particular the separation event, utmost priority was given to placing the Centaur close 
to its target, and sometimes this was at the expense of the SSC target accuracy. The 
separation event AV was used as a final maneuver to correct the Centaur trajectory 
in order to hit the target. The separation event was originally planned to be along the 
velocity vector but was taken off that in order to do a small correction thus causing a 
shift in the SSC trajectory. 

The largest source of error in the SSC impact location should have been due to the 
9 m/s braking burn that took place 40 minutes after separation. Other error sources 
include the large amount of attitude control thrusting during the various slews, as well 
as the tight dead band and more frequent quaternion updates used during observations. 
We estimated the braking burn error to be about -1% based on the post-burn OD (i.e., 
a slightly cold burn). The trajectory and image only analysis methods have small errors 
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Fig. 6 NSPl boresight angle from Centaur impact position in degrees as a function of time 

associated with their location estimates relative to the distance from the target and 
are in close agreement. The errors on the SSC trajectory as calculated by the JPL 
Orbital Determination team are 3 m x 75 m, Icr, where the 3m is in the Earth vector 
(as projected onto the lunar surface) and the 75 m is orthogonal to that. These are 
dominated by the errors of the tracking process. 



6 Pointing Accuracy 

The Near Infrared Spectrometer 1 (NSPl) was the primary water-detection instrument 
on LCROSS. Knowing the Centaur impact location allows the accuracy with which the 
SSC pointed this spectrometer at the impact to be evaluated. Figure 6 shows the angle 
between the NSPl boresight and the Centaur impact location over the last 15 minutes 
of the mission. The SSC kept the Centaur impact location continuously within the 0.1 ° 
of the NSPl field-of-view (FOV) from 5 minutes prior to Centaur impact until 3 minutes 
after. The angle remains under 0.5° until 11:35:07 UTC, (27 seconds before loss-of- 
signal). The SSC kept the Centaur impact location within the MIRl Camera fields-of- 
view until 11:35:32, (2 seconds before loss-of-signal) . The maximum error before the 
SSC starts to deviate from the Centaur impact location is 0.198°. 



7 Impact Angle 

There are two components to finding the impact angle of the Centaur with respect to 
the local terrain slope: (1) the velocity vector of the Centaur at impact, as measured 
with respect to the surface of a perfect lunar sphere; and (2) the slope at the impact 
point given the local terrain. 

The most accurate trajectory model that we have post impact was used to generate 
the local Centaur velocity vector at impact, using STK. At impact, the local lateral 
velocity components are -0.00217 °/s (= -0.06574 km/s) in latitude and 0.058947 °/s 
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(=0.165944 km/s) in longitude. Given the total velocity of 2.506885 km/s (or local 
altitude rate of 2.50052 km/s) at impact, this implies a latitude impact angle of - 
1.50271° and longitude impact angle of 3.79549°. This is obtained by projecting the 
trajectory onto a plane of constant longitude and measuring the angle between normal 
and trajectory in that plane. A negative angle is defined as the Centaur velocity being in 
the negative latitude direction. The total impact angle is thus 4.08295 ° from the normal 
to the lunar surface (or, 85.917° from horizontal) if the local slope was zero. Thus the 
Centaur was heading in a South-Easterly (mainly Easterly) direction at impact. 

LOLA topographical data were used to generate a local slope at the impact site 
by measuring the altitude at the impact site as well as at adjacent nearest neighbour 
points in the north-south axis and the east-west axis. The LOLA slope map is shown 
in Figure 7. 

The closest slope measurement (at 172 m from the impact location) had a slope 
of -0.14° in the East- West axis and 1.83° in the North-South axis, where a positive 
slope (x, y) dips (eastward, northward) - thus the total slope was 1.83°. The average 
of the 6 closest points (see Appendix D) were -1.00° in the East- West axis and 1.04° 
in the North-South, with the standard deviation being 2.20° and 0.80° respectively. 
Thus at the Centaur impact site the slope is most likely dipped towards the North and 




Fig. 7 LOLA Slope Map. The yellow cross indicates the Centaur Impact Location. The scale 
is linear with a slope of translating to ° and a slope of 0.6 to 30.96 °. 
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West, (mainly North), with slopes — 0.14± 2.2 and 1.83 ±0.8°. This does not account 
for any gradients with length scales < 400m. 

Combining the two measurements, and assuming a flat terrain at a < 400m scale, 
we have that the impact angle with respect to the local slope 0.33 it 0.8 ° in latitude 
and 3.66 =b 2.2° in longitude. 

In both directions the slope tends to cancel the local impact velocity. The overall 
surface impact angle is thus 3.67 =b 2.3 ° in a mainly East direction. Thus one possible 
explanation for why the images of the Centaur impact crater seen in the NIR2 camera 
show a greater amount of eject a in the Easterly direction that the impact angle resulted 
in an asymmetrical plume with more ejecta in that direction. On the other hand, this 
does not seem to be consistent with models and ground based experiments which show 
little variation of ejecta symmetry with impact angle if one used high impact angles 
(close to normal). In general, it is perhaps worth noting that this impact angle is 
somewhat steeper than the pre-encounter estimations used in impact models and lab 



et al. 



models (of order 70 degrees), see e.g. V. V. Shuvalov et al. (2008); Korycansky, D.G 



(2009), 



7.1 Impact Angle Error Analysis 

The errors are dominated by the local slope errors. Given that the closest measured 
slope is some 172 m from the Centaur impact location, and given that the slope varies 
by of order 2 ° based on translations of that scale between different LOLA slope mea- 
surements, in the region close to the impact point, the error margins on this basis are 
approximately 2°. The standard deviation of the 6 local slopes gives 0.8° in the N-S 
direction and 2.2 ° in E-W direction. Meanwhile the la errors in the angle deriving 
from the trajectory errors are negligible at < 0.01 °. 



8 Mission Requirements &z Lessons Learned 

Table [5] hsts six LCROSS mission requirements related to the Centaur and the SSC 
impact locations, the impact angles and observing geometry, and provides a summary 
of how well they were met based on analysis herein. 

Regarding PR J 4.2.3, the project's informal targeting goal was 1.75 km, which was 
also met with a wide margin. The approximately lOOm-level targeting accuracy of the 
mission may be of interest not only for future lunar purposes but also for Near-Earth 
Object (NEO) missions, some of which are of order this size. There is no obvious reason 
why a repeat mission with the same capability avionics and tracking systems would 
not be able to achieve a comparable targeting accuracy. 

With respect to PRJ4.1.7 and PRJ6.4.14-16, the SSC kept the Centaur impact 
location continuously within the 0.1 ° of the NSPl field-of-view (FOV) from 5 minutes 
prior to Centaur impact until 3 minutes after. Note that the commands controlling 
pointing were generated from an impact site prediction made several hours before the 
impact and the SSC stayed within 0.1 ° of that prediction. 

The analysis herein might suggest potential ways one may have improved the mis- 
sion design: 
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1. Centaur tracking improvements: one could have added some device such as a cor- 
ner cube to enable laser ranging, or a transponder, to have enabled more precise 
tracking of the Centaur post separation from the SSC. 

2. Camera pointing and distortion: several analyses could have benefitted from better 
ground based calibration using reference targets. 

3. Reaction wheels: one of the main sources of error on the spacecraft trajectory was 
that due to errors in the use of attitude thrusters, which could have been nulled 
through use of reaction wheels for attitude control. 

Naturally, these steps would most likely have had a cost and schedule implication 
which may have been prohibitive in this cost- and schedule-capped mission. Reaction 
wheels would have reduced dispersions and also some operational risks, but would have 
forced re-allocation of resources used elsewhere, increasing risks there, etc. 

It can equally be stated, that each of the mission objectives, and in particular the 
central scientific questions, were answered with the mission as was. In fact one can ask 
whether LCROSS was over-designed in some areas and speculate on changes to the 
mission that may have reduced complexity and cost: 

1. Avionics: potentially a simpler/less capable avionics system could have been em- 
ployed since the targetting was met with such a high margin 

2. Tracking accuracy: the tracking accuracy could have been less and still met the 
mission requirements, (e.g. DSN could have been used less). As specific examples, 
mission operations could have foregone the update to the pointing quaternions that 
was made ^ 5 hrs before impact and the mission might have been able to use an 
earlier OD solution to generate the pre-separation load of quaternion and braking 
burn parameters in order to give more time to generate and do quality assurance 
on products. 



Table 5 Six key related mission requirements of LCROSS how they were met 



Req. No. Requirement Text 



Result 



LPRP4.4 The LCROSS impactor shall impact the 
target location at an angle greater than 60 
degrees to the plane of the lunar surface. 

PRJ4.2.3 The LCROSS shall be able to target a 10 
km radius impact area (3-sigma). 

PRJ6.4.8 LCROSS shall be capable of observing im- 
pact target location and the sunlit ejecta 
cloud with at least 0.1 Hz frequency from 
EDUS impact minus 5 minutes until EDUS 
impact minus 3 minutes 

PRJ6.4.14 The FS shall maintain payload pointing to 
the Centaur impact site at the time of Cen- 
taur impact with an accuracy of +/- 0.5 
degrees 

PRJ6.4.15 The FS shall maintain payload pointing to 
the Centaur impact site at Centaur impact 
plus 60 seconds with an accuracy of +/- 
0.64 degrees. 

PRJ6.4.16 The FS shall point the nadir viewing NIR 
Spectrometer to the Centaur impact site at 
a time of 60 seconds prior to S-SC impact 
to an accuracy of +/- 5 degrees. 



Impact angle 85.92 ± 0.01° 
w.r.t. local slope. 

Distance to target: 146 m, 
max. mean plus 3<t range < 
500m 

Sample rates varied by instru- 
ment at time, but all were < 
0.1 Hz over this period and 
were observing the target loca- 
tion and ejecta cloud. 
Pointing error < 0.1° at im- 
pact 



Payload pointing to impact lo- 
cation < 0.1° at impact plus 
60s 

Payload pointing to impact lo- 
cation < 0.2° at SSC impact 
-60 s 
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3. Fuel for ACS: 

(a) less fuel could have been used for ACS while still meeting the pointing require- 
ments by allowing larger attitude control margins of error 

(b) less fuel could be used by not doing force balancing while still meeting the 
spacecraft targeting requirements 

4. Fuel for TCMs: perhaps fewer Trajectory Correction Maneuvers (TCMs) could have 
been made 



9 Conclusion 

We hope that this paper supports understanding the LCROSS science data. In partic- 
ular that: 

1. The impact imagery registration and locations could aid follow-up observations 
from the ground and from lunar orbit, particularly from LRO e.g. for example 
the LRO-LOLA instrument to attempt more careful study of the topography sur- 
rounding the impact locations to try to identify the crater (and thus independently 
measure features); 

2. The impact location, derived from several independent methods, may help to con- 
strain or improve the map tie errors in the lunar coordinate reference systems. 

3. Full registration of all descent images from the automated bundle adjustment pro- 
cess used herein could provide additional, useful data to users of the LCROSS data 
and will be provided to the PDS; 

4. The impact angle of the Centaur with respect to the terrain may aid the develop- 
ment of impact cratering and ejecta process models; and 

5. The pointing accuracy and motion of the LCROSS spectrometer boresight may 
help users of the data to better understand the spectra. 

LCROSS mission requirements relating to targeting, instrument fields of view, 
pointing accuracy and impact angle were all met. In fact 5 of the 6 did so by bet- 
ter than a factor of 5. 



10 Acknowledgements 

We would like to acknowledge the efforts and support of the following teams/entities 
without whom this analysis would not have been possible: the Mission Operations Sup- 
port team at NASA- Ames Research Center (ARC), the Mission Design Team at ARC, 
the Maneuver Design Team at NASA-Goddard Spaceflight Center (GSFC), the Orbit 
Determination Team at the Jet Propulsion Laboratory, the Subsystem analysts and 
engineers at Northrup Grumman Space Technology and Northrup Grumman Technical 
Services and the operators of the DSN facilities at Goldstone, Madrid and Canberra. 



References 

Andrews, D., LCROSS Lunar Impactor - Lessons Learned from a Small Satellite Mission, 61st 
International Astronautical Congress, Prague, CZ, IAC-10-B4.8.2 (2010). 



16 



Arganda-Carreras, I., Sorzano, C. O. S., Marabini, R., Carazo, J. M., de Solorzano, C. O., 
and Kybic, J., Consistent and Elastic Registration of Histological Sections using Vector- 
Spline Regularization, Lecture Notes in Computer Science, CVAMIA: Computer Vision 
Approaches to Medical Image Analysis, (Springer Berlin / Heidelberg) 4241, 85 (2006). 

Broxton, M., Nefian, A. V., Moratto, Z., Kim, T., Lundy, M., Segal, A. V. 3D Lunar Terrain 
Reconstruction from Apollo Images, Proceedings of the 5th International Symposium on 
Advances in Visual Computing: Part /, 710 (2009). 

Colaprete, A., Schultz, P., Heldmann, J., Wooden, D., Shirley, M., Ennico, K., Hermalyn, B., 
Marshah, W., Ricco, A., Elphic, R. C, Goldstein, D., Summy, D., Bart, G. W., Asphaug, 
E., Korycansky, D., Landis, D., Sollitt, L., Detection of Water in the LCROSS Ejecta 
Plume Science, 330 (6003), 463 (2010) 

Cooley, D.S., et al.. Mission Design for the Lunar CRater Observation and Sensing Satellite 
(LCROSS), AIAA/AAS Astrodynamics Specialist Conference, Toronto, Canada, (2010) 
(to appear). 

Ennico, K., et al., LCROSS Science Payload Ground Development, Test and Calibration Re- 
sults, 39th Lunar and Planetary Science Conference, Lunar and Planetary Science Insti- 
tute, League City, TX, (1391), 1474 (2008). 

Ennico, K., Colaprete, A., Shirley, M. and D. Wooden, Lunar Crater Observation and Sensing 
Satellite (LCROSS) Instrument Calibration Summary, jfill in PDS reference^ (2010). 

Ennico, K., M. Shirley, Colaprete, A. and Osetinsky, L., The Lunar Crater Observation and 
Sensing Satellite (LCROSS) Payload Development and Performance in Flight, In This 
Issue. 

Korycansky, D. G., Plesko, C.S., Jutzi, M., Asphaug, E., and Colaprete, A.. Predictions for 
the LCROSS mission Meteoritics and Planetary Science, 44 (4), 603 (2009) 

Heldmann, J., Marshall, W., Colaprete, A., Goldstein, D., and Summy, D., Evolution of the 
dust and water ice plume components as observed by the LCROSS Visible Camera and 
UV-Visible Spectrometer [In preparation]. 

LCROSS Flight Director's Blog, ,http://blogs.nasa.gov/cm/newui/blog/| view- 
postlist.jsp?blogname= Icrossfdblog (2009^ 

Moratto, Z., et. al, LCROSS Derived Map Products, 3rd Annual Lunar Science Forum, NASA- 
Ames, USA, Poster 220, (2010) 

Shuvalov, V.V. and Trubetskaya, LA., Numerical Simulation of the LCROSS Impact Experi- 
ment Solar System Research, 42 (1), 1 (2008) 

Schultz, P., Hermalyn, B., Colaprete, A., Ennico, K., Shirley, M., and Marshall, W. S., The 
LCROSS Cratering Experiment Science, 330 (6003), 468 (2010) 

Strong, J., et al.. Transport and Use of a Centaur Second Stage in Space, AIAA Space Ops 
2010 Conference, Delivering on the Dream, AIAA, Washington, DC, Paper 2010-2197 
(2010). 

Tompkins, P., et al.. Flight Operations for the LCROSS Lunar Impactor Mission, AIAA Space 
Ops 2010 Conference, Delivering on the Dream, AIAA, Washington, DC, Paper 2010-1986 
(2010). 

Tompkins, P., et al.. Flight Team Development in Support of LCROSS: A Class D Mission, 
AIAA Space Ops 2010 Conference, Delivering on the Dream, AIAA, Washington, DC, 
Paper 2010-2223 (2010). 



17 



A Appendix A: Source Data and Analysis Process 

The data sources used for all analyses herein are described in overview in Table [6] The process 
flow, showing input data used for each impact location estimation process shown schematically 
in Figures 8 and 9 for the SSC and Centaur respectively. 



Table 6 Source data for location measurements. Note that the LOLA data was produced in 
October 2011 based on > 20 months of integrated data and that the resolution of the DSN is 
Im radial to the antenna (in the Earth direction) and 75m in orthogonal directions 



Type 



Source Description 



Area Resolution Time Stamp 

(km) (m) 



Terrain 


LRO 


LOLA terrain maps 


10 X 10 


~ 40 


Jan. 2010 


NIR 


LCROSS 


NIR image 


4x5 


~ 10 


11:35:31.404 


MIR 


LCROSS 


MIR image 


2x3 


~ 20 


11:35:31.401 


Radar 


Goldstone 


Radar images before 


5x7 


~ 5 


May & Nov. 






&; after 






2009 


Ranging/ 


DSN 


SSC ranging data 


n/a 


1 X 75 


~ 1/hour until 


Doppler 




taken up until l.Shrs 




10:17:00.000 






before impact 








Attitude 


LCROSS 


SSC attitude sensors 


n/a 


< 0.1° 


~ IHz through 






data via telemetry 






descent 



/ SSC 
Ranging y 



SSC 




Trajectory 





SSC 
/ Tra^ctoryy 



Trajectory only 



Fig. 8 Flow diagram of SSC impact determination process. On the LHS are the input data. 
On the RHS are the impact coordinate products. In the middle are the intermediate processes 



Image only 




Hybrid Methods 
Trajectory only 



Fig. 9 Flow diagram of Centaur impact determination process. On the LHS are the input 
data. On the RHS are the impact coordinate products - the four semi-independent method- 
ologies (the image only and trajectory only are completely independent, and the manual and 
automated hybrid methods are semi independent of the other two). In the middle are the 
intermediate processes. 
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B Appendix B: Further Details on The Centaur Impact Location 

In this section, each of the four Centaur impact location determination methods, as summarised 
in Table 2, are detailed. 



B.l Centaur Trajectory Only 

LCROSS tracking and orbit determination was performed by JPL using transmitter ranging 
and Doppler data from the Deep Space Network (DSN). With no onboard transmitter, the 
Centaur was untracked from the time it separated from the SSC until impact 9.7 hours later. 
Any trajectory-based estimate of the impact location must address the trajectory perturbation 
introduced by the separation event. Two methods of estimating the separation perturbation 
were used yielding two estimates of the post-separation Centaur state vector. These state 
vectors were then propagated forward until they intersected the surface. 

The first method involves using the results of pre-separation orbit determinations to calcu- 
late the state vector of the combined SSC/Centaur at the time of separation, then perturbing 
that vector by the impulse applied (to each object) by the separation springs, as known from 
ground calibration, to derive the Centaur's post separation state vector. This estimate of the 
separation impulse is purely analytical. 

The last state for the SSC prior to separation was calculated at 09 Oct 2009 00:20:00 UTC. 
The Cartesian coordinates in J2000 with the Earth as the central body were: 

Position [km]: 87774.263494 340013.79800, 122539.97427 
Velocity [km/s]: -0.84355753755, -0.09197672465, 0.85811369196 

The separation event AV target was 0.721 m/s with a 3a uncertainties of 9 mm/s, per axis. 
The target attitude quaternion (in J2000) of the spacecraft at separation was: 0.36790183411, 
0.73856335829, -0.40935830939, 0.38935611048. The spacecraft attitude quaternion as mea- 
sured by the on-board attitude determination system at the time closest to separation (01:49:58.67) 
was -0.368135433852522, -0.738446417838845, 0.40900543141651, -0.389727786275225, suggest- 
ing an error of < 0.066 ° from the command (target) quaternion. The resulting Centaur post- 
separation state vector at 09 Oct 2009 01:50:00 is: 

Position [km]: 83216.293080, 339463.90550, 127193.59414 

Velocity [km/s]: -0.84438148975, -0.11192890266, 0.86605975234 

The second method involves using tracking data from before and after separation plus 
Newton's Third Law. Here the trajectory of the combined SSC-Centaur vehicle prior to sep- 
aration, together with the trajectory of SSC post-separation (but prior to braking burn) are 
taken, along with the SSC and Centaur relative masses at separation. Together these enable 
one to deduce the post-separation Centaur trajectory based on its equal and opposite reaction 
at the separation event. The first state for the SSC post-separation was calculated for the time 
09 Oct 2009 02:28:59. The Cartesian coordinates in J2000 with the Earth as the central body 
were: 

Position [km]: 81240.391374, 339191.93319, 129222.82223 
Velocity [km/s]: -0.84480650128, -0.12083383916, 0.86981590338 

The SSC mass at separation was 617.18 kg and the Centaur mass at separation 2271.61 
kg. This, together with the SSC-Centaur pre-separation state already given above for the first 
method, were propagated to the separation time of 09 Oct 2009 01:50:00 - in the case of the 
SSC post-separation state the propagation being backwards in time. The resulting Centaur 
post separation state vector, assuming an instantaneous maneuvre, and at the same time is: 

Position [km]: 83216.275714, 339463.90071, 127193.61830 
Velocity [km/s]: -0.84444756135, -0.1118950345, 0.86592455048 

The two post-separation Centaur state vectors were propagated forward until they in- 
tersected the lunar surface using STK-Astrogator (9.0). This propagator included a 70x70 
LP-based lunar gravity model, the Sun and Earth as point masses, as well as solar radiation 
pressure. Based on LOLA data, the surface altitude of -3.82633 km from a mean lunar radius 
of 1737.4 km. The JPL-based Orbit Determination team replicated the second method using 
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Table 7 Summary results of the trajectory-only impact position estimates. Altitude is with 
respect to a sphere of radius 1737.4 km. Latitude and longitude are given in Mean Earth, 
selenocentric coordinates. 



Source 


Latitude (°) 


Longitude ( ° ) 


Altitude (km) 


Trajectory 1 


-84.6749921 


-48.7001 


-3.82817 


Trajectory 2 


-84.675 


-48.703 


-3.82633 


Centaur JPL 


-84.6749842 


-48.6972 


-3.83 


Average 


-84.6750 


-48.7003 


-3.82862 


Average + outgassing 


-84.6827 


-48.6688 


-3.82862 



an alternative propagation program to produce a third, independent trajectory-based impact 
location estimate. The predicted impact locations from these methods and their average are 
below. 

These estimates do not include a potentially significant, intermittent perturbation: out- 
gassing from ice remaining on the dark side of the Centaur. Deposited by the moist Florida 
air before launch, trajectory perturbations from ice outgassing were observed during a ma- 
neuver in August, two months after launch and one month before the impact. In the nominal 
LCROSS attitude, this ice is shadowed by the Centaur and inert. But after separation, the 
Centaur tumbled, exposing the ice to the sun and causing off gassing and consequently a small 
but unmeasured perturbation. 

Dispersions due to ice outgassing will not be perfectly circular. Since the outgassing occurs 
only when the ice is lit, it tends to push the trajectory anti-sunward. This effect is more 
pronounced when the Centaur rotation period is less than the time before impact, as it was. 
Therefore, we expect outgassing to cause the Centaur to land anti-sunward of the trajectory- 
based prediction point. 

The outgassing acceleration observed in August was 5.7 mm/s per hour, applied to the 
mated SSC/Centaur (~ 3150 kg). Acceleration of the Centaur alone (2366 kg) would thus be 
approximately 7.6 mm/s per hour. If these perturbations are assumed to be applicable for half 
of the post-separation time (9.7 hours), then the total deflection is 640 meters. After applying 
cosine losses to accelerations during Centaur rotation, the impact deflection can be bounded to 
be > 410 m. The lower limit is 100m. Thus the average would be 250m with an approximately 
50m 1(7 error to account for these 3a bounds. These are the dominant contributions to the 
errors given in the main text (the propagation and initial condition errors being smaller). The 
results from these methods are summarised in Table 7. 



B.1.1 Centaur Trajectory Only Error Analysis 

The standard deviation of the trajectory methods is 89m in latitude and 107m in longitude. 
However, the standard deviation of this clustering does not include any systematic offset of 
this methodology. To understand the true error, one needs a bottom up approach to error 
estimation. 

The determination and prediction of the LCROSS trajectory was performed by the Orbit 
Determination (OD) team at JPL. Based on a solution containing tracking data up to the sep- 
aration event, JPL's assessment of the Icr knowledge of the Centaur Impact is a nearly circular 
ellipse with radii of 94 meters. This is primarily attributable to three sources of uncertainty: (1) 
the state knowledge at the time of the last tracking data observation of the Centaur trajectory 
(i.e. just before separation). Speciflcally the cross-track knowledge at 36 meters and 2.3 mm/s, 
1(7, equivalent to a deflection of 88m; (2) the effects of the separation AV uncertainties on 
the relatively massive Centaur. Speciflcally these were 3 mm/s, per axis, Icr, equivalent to a 
deflection of 105m; and (3) the errors of potential outgassing as residual water ice is burned off 
whenever the Centaur dark side tumbles into the sunlight. The uncertainties on this amount 
to 50m 1(7. Note that these need to be applied around a position that is offset by the mean 
expected offgassing. Adding these errors together we flnd a total error of 141 x 141m. 
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B.2 Centaur Imagery Only 

The method for locating the impact based on imagery only uses the co-registration of impact 
crater images described in Sectionjs] Since LOLA data is registered (in the Mean Earth frame), 
the process of registering the NIR2, MIRl and GSSR images (which show the Centaur crater) 
to LOLA allows us to estimate the Centaur impact location in absolute terms. The crater 
feature is most clearly visible in the NIR2 image. The location of the centre of this feature is 
-84.6774°, -48.691 ° in Mean Earth, selenocentric coordinates. 

Note that with the data used herein, it has not been possible to identify the Centaur crater 
in the LOLA data but further data may allow the crater to be resolved. Note also that the 
feature in the NIR2 imagery which is over 150 m in diameter, is dominated by the bright ring 
of the freshly ejected material around the Crater, which has different scattering properties, 
but is likely to make little difference topologically, so the LOLA feature would be expected to 
be considerably smaller, given the crater size of 25-30m ( |P. Schultz et al.|2010| ). 

B.2.1 Image Only Method Error Analysis 

The error for this approach is given by three factors. First the error in the image registration. 
This was estimated by repeating the whole registration process and measuring the difference 
in impact location. The maximum range of trials was found to be ~ 100m in latitude and 
longitude. Second the error associated with measuring the impact site in the NIR2 image. 
This was estimated by taking repeated measurements of a given registered image set, outlining 
the errors. Further, depending on where the impact point is within the features, the errors 
build. The NIR2 feature is 100x120m. The overall error associated with where the impact 
point in these features is approximately 70x70m. Finally, there is the accuracy with which the 
LOLA data has been registered to the lunar coordinate frame (the map tie error), which is 
approximately 40m. These errors are independent so we combine them using a quadratic mean 
to give 130m. This is approximately equal in both axes. 



B.3 Hotspot Forward Projection (Manual) 

Two more approaches were developed which used a combination of imagery and tracking data 
to determine the impact location. The first involves using Satellite ToolKit (STK) to project 
the vector through the hotspot in the MIRl images from the SSC (knowing it's attitude and 
position) to the surface. The second method is a semi-automated image registration process 
that was applied to all of the images. This process was also used to generate image registration 
data to be contained in supplemental data to the LCROSS delivery to the Planetary Data 
System (PDS). 

The Centaur impact location is clearly visible in MIRl (thermal) images taken in the last 
30 seconds of the mission. For each image, the first technique projected a ray from the SSC 
centre's position when the image was captured, through the hotspot in the image plane of 
the camera to the surface. This produced a cluster of points (corresponding to different MIRl 
image times) on the surface whose locations were averaged to estimate the Centaur impact 
location. The information needed to implement this method is: (1) the SSC trajectory; (2) 
the SSC attitude history; (3) the MIRl camera attitude relative to the SSC body coordinate 
system; (4) the pixel locations of the hotspots in the MIRl images; and (5) the altitude of the 
lunar surface terrain near the impact location. 

This approach was implemented using STK and the LCROSS DataBrowser. STK had 
already been used to model the SSC trajectory and attitude. Within STK, the Moon was 
modelled as a sphere of radius 1737.4 (lunar mean radius) with a texture generated from 
LOLA terrain data. Also, for the purposes of this analysis, the area surrounding the Centaur 
impact was modelled as being flat at an altitude of -3.82633 km from the lunar mean radius. 
Therefore, STK took care of items 1, 2 and 5. 

The LCROSS DataBrowser is custom software that we use to explore payload imagery and 
spectra. Here, it's simply a program for selecting between MIRl images along a timeline and 
displaying them in a window. The DataBrowser was modified so that the window displaying the 
images could be made partially transparent so that STK windows could be seen underneath. 
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Table 8 Impact location estimates resulting from one iteration of the manual image registra- 
tion method. Latitude and longitude are given in Mean Earth, selenocentric coordinates. 



MIRl image 


Latitude (°) 


Longitude ( °) 


Measurement 1 


-84.6820039 


-48.74303581 


Measurement 2 


-84.68174335 


-48.74546489 


Measurement 3 


-84.68177716 


-48.75053408 


Measurement 4 


-84.68163383 


-48.75641848 


Measurement 5 


-84.68171144 


-48.75977728 


Average 


-84.68177394 ± 0.000139 


-48.75104611 ±0.00708 



Because no images show the Centaur crater together with enough context to match against 
the LOLA-derived image, the camera pointing was generated from an image taken ~ 20 minutes 
earlier and then reused for later images that show the Centaur crater hotspot. The pointing 
was generated by selecting an image in the DataBrowser, transferring the time to STK to 
generate an image from the S-S/C position and attitude at that time, then adjusting the both 
the STK image (varying rotation around the boresight) and the DataBrowser image (adjusting 
X and Y scale and position) so the images matched. 

This process generates the pointing for the camera implicitly, i.e., it doesn't generate the 
pointing as a quaternion or transformation matrix, but it does align the two programs so that 
images at other times will remain aligned. Once this has been done, later images showing the 
hotspot can be superimposed over the STK display. Within STK, a target object was created 
at an altitude of -3.82633 km and moved in latitude and longitude to lie under the hotspot. 
This process was repeated for each of five MIRl images resulting in 5 targets whose positions 
were averaged. The entire process was then repeated from the beginning. The results from 
these methods are summarised in Table 8. 



B.3.1 Hotspot Forward Projection Error Analysis 

Errors in the Hotspot Forward Projection Method arise from (1) the SSC trajectory, (2) the 
SSC attitude, (3) MIRl camera pointing, and (4) alignment of surface targets with the hotspot 
in the MIRl images. Errors in camera pointing arise from (5) manual image alignment and (6) 
registration of the LOLA terrain, to which an MIRl image was aligned, to STK's Mean Earth 
coordinate frame. #6 was trivial compared to #5 because the SSC was at a high altitude when 
the MIRl image was taken. 

We divide these into errors arising in the trajectory estimate and all errors arising in 
the image alignment and attitude. The different image-related error sources are complex, so 
we estimate the error by assuming a 2-D normal distribution for the resulting impact crater 
locations. From the 10 samples, the Icr error is 18 m. The SSC trajectory error is 3 m x 75 m, Icr 
in the lunar surface plane, where 3 m is in the Earth vector, and the other is orthogonal. These 
two error sources are independent, so we combine them with their quadratic mean (RMS), 
handling each axis separately. This yields an error ellipse for manual image registration of 21 
m X 77 m. Including the 40 m LOLA map-tie error, this goes to 45 m x 87 m, Icr, where the 
smaller error is along the Earth-axis and the larger is orthogonal to that. 



B.4 Image Registration by Bundle Adjustment (Automated) 

The Bundle Adjustment (BA) approach was applied to solve for the relative pose and pro- 
jective properties of the five LCROSS cameras with respect to the SSCs frame. With that 
solution, all LCROSS images where projected onto the LOLA topography where the latitude 
and longitu de of individual pixe ls could be extracted. Implementation and algorithm were 
taken from jBroxton et al.|2010| >. 

This algorithm reduces the projection error of 3D features into cameras by comparing the 
observed location of features in the images compared to where the camera models predicts 
them. The location of the 3D features were solved simultaneously along with the the rotation 
of the cameras with respect to SSC and their individual focal lengths. Two types of features 
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where used in the optimization, Tie Points and Ground Control Points. Tie points are features 
found in 2 or more images where only the their pixel location was recorded. Their 3D location 
was seeded with a a mid-point triangulation method. Ground Control Points (GCPs) are image 
measurements similar to Tie Points that also have been registered to LOLA topography by 
hand. GCPs allow for an additional cost function in BA where their solved 3D location must 
stay close to their observed location on LOLA. 

Two time instances where found where the five cameras (VIS, NIRl, NIR2, MIRl and 
MIR2) were all taking exposures within a couple of seconds. 9 GCPs were tracked in those 10 
images simultaneousl y. An addition al 58 Tie Points were tracked but with varying coverage 
among the 10 images ( |Moratto|2010| >. The results from these methods are summarised in Table 
9. 



Table 9 Summary of measurements of impact location by automated bundle adjustment. 
Latitude and longitude are given in Mean Earth, selenocentric coordinates. 



MIRl image 


Latitude (°) 


Longitude ( °) 


1 


-84.676776 


-48.7064 


2 


-84.676755 


-48.7036 


3 


-84.67627 


-48.7098 


4 


-84.677526 


-48.7001 


5 


-84.676863 


-48.7177 


6 


-84.676997 


-48.7006 


7 


-84.677283 


-48.7031 


8 


-84.676841 


-48.7149 


9 


-84.680077 


-48.7204 


10 


-84.679855 


-48.7263 


11 


-84.679866 


-48.7396 


12 


-84.679064 


-48.742 


13 


-84.680799 


-48.7561 


14 


-84.680985 


-48.7594 


Average 


-84.6783 ±0.0017 


-48.721 ±0.020 


la error (m) 


52 


57 



B.4-1 Error Analysis 

Errors and uncertainty where taken from derivatives of the last iteration of Bundle Adjustment. 
These translate to a ground position error of 52 x 57 m when projecting images near impact. 
In addition, LOLA la map tie errors are approximately 40m. The trajectory errors are 3x75m. 
Thus a total estimate might be 66xl02m Icr, where the smaller axi s is in th e Earth vector 
(projected onto the lunar surface) and the larger is orthogonal (Moratto|2010| 
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C Appendix C: Further Details on the SSC Impact Location 

In this section, the SSC impact location determination method based on the trajectory only 
is detailed. The SSC was tracked until quite close to impact, providing a fairly accurate loca- 
tion. The last state vector available for the SSC prior to impact is for 09 Oct 2009 10:17:00 
(approximately 79 minutes before impact). The Cartesian coordinates in J2000 with the Earth 
as the central body are: 

Position [km]: 57759.363734, 333782.02444, 155582.61612 
Velocity [km/s]: 0.798778319, -0.327045732, 1.16015291 

This state was propagated by STK-Astrogator (9.0) until the trajectory met the lunar 
surface. Astrogator used a 70x70 LP-based lunar gravity model, the Sun and Earth as point 
masses, and included solar radiation pressure and general relativity. The altitude cut off was 
taken to be -3.82633 km from a mean lunar radius of 1737.4km, based on LOLA data. The 
impact position determined for the spacecraft was: Latitude, Longitude, altitude: -84.719 °, 
49.610 °, -3.80909 km. 

The errors on the SSC trajectory as calculated by the JPL Orbital Determination team 
are 3 m x 75 m, la, where the 3m is in the Earth vector (as projected onto the lunar surface) 
and the 75 m is orthogonal to that. These are dominated by the errors of the tracking process. 
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D Appendix D: Local Slope Measurements 

Table 10 shows slope measurements made by LOLA that are local to the Centaur impact site. 
These are used in Section 7. 



Table 10 Local Slope measurements. (Note the impact location is -84.6775, -48.7216) 



Slope 


Latitude (°) 


Longitude ( ° ) 


E-W slope 

n 


N-S slope 

n 


Impact Dis- 
tance (km) 


Point 1 


-84.6566 


-48.75339 


-3.01871 


0.691863 


0.649182 


Point 2 


-84.6615 


-48.75 


-0.01793 


-0.24683 


0.497813 


Point 3 


-84.6676 


-48.76142 


-4.26773 


1.646793 


0.321277 


Point 4 


-84.6726 


-48.75794 


-0.14788 


1.828963 


0.172431 


Point 5 


-84.69 


-48.77573 


-0.2415 


1.588041 


0.38113 


Point 6 


-84.6944 


-48.76999 


1.67379 


0.71209 


0.50838 


Average 

Standard 

Dev. 


-84.6738 
0.015325 


-48.76141 
0.0098531 


-1.00333 
2.198553 


1.03682 
0.796945 





